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ABSTRACT 



Aims. We analyse the oscillatory properties of resonantly damped transverse kink oscillations in two-dimensional prominence threads. 
Methods. The fine structures are modelled as cylindrically symmetric magnetic flux tubes with a dense central part with prominence 
plasma properties and an evacuated part, both surrounded by coronal plasma. The equilibrium density is allowed to vary non-uniformly 
in both the transverse and the longitudinal directions. We examine the influence of longitudinal density structuring on periods, damping 
times, and damping rates for transverse kink modes computed by numerically solving the linear resistive magnetohydrodynamic 
(MHD) equations. 

Results. The relevant parameters are the length of the thread and the density in the evacuated part of the tube, two quantities that 
are difficult to directly estimate from observations. We find that both of them strongly influence the oscillatory periods and damping 
times, and to a lesser extent the damping ratios. The analysis of the spatial distribution of perturbations and of the energy flux into the 
resonances allows us to explain the obtained damping times. 

Conclusions. Implications for prominence seismology, the physics of resonantly damped kink modes in two-dimensional magnetic 
flux tubes, and the heating of prominence plasmas are discussed. 
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1. Introduction 

Quiescent filaments/prominences are cool and dense magnetic 
and plasma structures suspended against gravity by forces 
thought to be of magnetic origin. In spite of their physical prop- 
erties, with temperatures and densities that are akin to those 
in the chromosphere, some as yet not well determined mecha- 
nisms provide the required thermal isolation from the surround- 
ing coronal plasma and mechanical support during typical life- 
times from few days to weeks. The magnetic field that pervades 
these structures is believed to play a key role in the nature and 
the thermodynamic and mechanical stability of prominences. 
Early observations carried out with good seeing c onditions 
pointed out that prominences consis t of fin e threads (de Jager, 
1959; Ku perus & Tandberg-Hanss en. 1967). More recent high- 
resolution H a observations obtai ned wi t h the Swedish Solar 
Telescope (SST) in La Palma dLin et all 120051) and the Dutch 
Open Telescope (DOT) in Tenerife (Heinzel & Anzer, 2006) 
have allowed to firmly establish the filament sub-structuring and 
the b asic geometrical and physical properties of th reads (see 
also lEngvoldlll998l:lLin et all I2005L 120081: iLinL 120 1 Dh . The sub- 
structure of quiescent prominences is often composed by a myr- 
iad of horizontal, dark and fine threads, made of cool absorbing 
material, believed to outline magnetic flux tubes (|Engvoldl|l998l 
2008: lLh3.l2004l:lLin et all 120051 12008: Ma rtin et alll2008h . The 
tubes are only partially filled with cool and dense plasma and 
their total length is probably much larger (~ 10 5 km) than the 



threads themselves. The measured average width of resolved 
threads is about 0.3 arcsec (~ 210 km) while their length is be- 
tween 5 and 40 arcsec (~ 3500 - 28 000 km). T he absorbing coo l 
material is usually visible for up to 20 minutes dLin et all 12005). 
The measured widths are close to the current resolution limit, 
~ 0.16 arcsec at the SST, hence thinner structures are likely to 
exist. 

Small amplitude oscillations in pr ominence threads a re fre- 
quently observed (see reviews by Oliver & Ballester, 2002; 
l Engyoldl 12004 IWiehi 12004 iBallestetl 120061: iBaneriee et all 
120071: lEngvoldl I2008t lOliverl l2009t iBallesterl |20 1 Ol) . Ear ly two- 
dimensi o nal o bservations of filaments ( Yi & Engvold, 1991; 



lYi et all 1199 lh revealed that individual threads or groups of 
them oscillate with periods that range between 3 and 20 min- 
utes. Recent relevant examples are traveling waves propagat- 
ing along a number of threads with average phase speed of 12 
km s , wavelength of 4 a rcsec, and oscil latory periods that 
vary from 3 to 9 minutes dLin et al ., 2007), the both propa- 
gating and stan ding oscillations dete cted over large areas of 
prominences by iTerradas et all d2002l) and iLinl d2004l) . as well 
as obse rvations from instruments onboard space-crafts, such as 
SoHO dBlanco et all 1 19991: iReenier et a 1 l200ll iPouget et al " 



2006) and Hinode (Oka motoet all 120071: 



Terradas et al., 2008b; 



mget 



Ning et al., 2009). The transverse oscillation nature of some of 
these events has been clearly established by iLin et a D d2009h 
by combining H a filtergrams in the plane of the sky with H a 
Dopplergrams which allow to detect oscillations in the line- 
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of-sight direction. A recurrently observed property of promi- 
nence oscillations is their rapid temporal damping, with per- 
turbations d ecaying in time-scales of only a few oscilla- 
tory peri ods (Landman et all Il977t iTsu baki & Takeuchi, 1986; 
Tsub akjLll988HWiehr et allll989tlMolownY-Horas et all fl999t 



Terradas et alll2002tlLinll2004tlNing et alll2009l) . 

Transverse thread oscillations are commonly interpreted in 
terms of standing or propagating magnetohydrodynamic (MHD) 
kink waves. The measured periods are of the order of a few 
minutes and t he wavelength s are i n between 3000 - 20 000 
km, although Oka moto et al.l d2007l) report larger wavelengths 
more consistent with the standing wave interpretation. The mea- 
sured wave quantities allow us to derive phase speeds that 
are consistent with the kink speed in magnetic and plasma 
configurations with typical properties of prominence plasmas. 
The MHD wave interpretation of thread oscillatio ns has al- 
lowed the development of the oretical models (se e Ballester, 
120051 1 20061 for recent reviews) . Hoarder et all dl997h : |Piaz etalj 



(2001 



2003) considered the MHD eigenmodes supported by a 



filament thread modelled in Cartesian geometry. More realis- 
tic studies using cyli ndrical configurat i ons have extended the 
initial investigations (iDfaz et all 120021 : iDvmova & Rudermanl 
120051: iDfaz & RobertsL 120061) . These studies have determined 
the frequencies and confinement properties of the perturba- 
tions as functions of the length and the width of the threads. 
Theoretical damping mechanisms have also been developed (see 
iBallesteri [20101: lArregui & Ballesteil 120101, for recent reviews). 
A systematic com parativ e study of different mechanisms has 
been presented by Solej] d2010h . who assesses the ability of 
each mechanism to reproduce the observed attenuation time- 
scales. The considered mechanisms include non-ideal effects, 
such as radiation and thermal conduction, partial ionisation 
through ion-neutral collisions, ion-electron collisions, and res- 
onant absorption due to coupling to Alfven and slow waves. 
Non-ideal effects do not seem to pr ovide t he required attenuation 
time-scales for kink oscillations dBallail [20 03; Carb onell et all 
2004: |Terradas et alll200W2005l:ISoler et al.Ll2007ll2008h . lSolerl 
(2010) finds that resonant absorption in the Alfven contin- 
uum is the only mechani sm able to produce the o bserved at- 
tenuation time-scales (see ISoler et all I2008L l2009allbllcl for de- 
tails). Resonant dam ping was first considered in this context by 
Arre gui et al.l (2008), and has been subsequently studied in com- 
bination with damping in th e slow continuum and partial ioni- 
sation bv lSoler et alj d2009alfbllcl) . These studies have considered 
one-dimensional models with a density variation in the trans- 
verse direction only. On the other hand, theoretical studies that 
take into account the longitudinal density structuring, hence in- 
corporating the fact that magnetic tubes supporting threads are 
only partially filled with cool plasma, consider piece-wise ho- 
mogeneous models in the transverse direction, and hence rule 
out resonant damping. 

A recent investigation by Soler et al.lJ20 10) is the exception. 
These authors have obtained analytical and semi-analytical ap- 
proximations for periods, damping times, and damping ratios of 
standing kink modes in a two-d imensional promi nence thread 
model. However, the analysis bv lSoler et al.l d2010h is restricted 
to the thin tube and thin boundary approximations, radial inho- 
mogeneity is constrained to the dense cool part of the tube only, 
and the density in the evacuated part of the tube is taken equal to 
the coronal density, further limiting the prominence threads that 
can be theoretically modelled. These limitations are removed 
in the present investigation. We adopt a fully inhomogeneous 
two-dimensional prominence thread model. The density distri- 
bution is allowed to vary non-uniformly in both the transverse 



and longitudinal directions. We also add another relevant param- 
eter with seismological implications, namely, the density in the 
evacuated part of the tube. By combining the two parameters that 
characterise the longitudinal structuring in prominence threads, 
the length of the thread and the density in the evacuated part of 
the tube, a wide range of prominence threads with very different 
physical conditions can be modelled and their oscillatory prop- 
erties characterised. In addition, a general parametric study is 
numerically performed, that allows us to go beyond the thin tube 
and thin boundary approximations. 

Besides obtaining the influence of longitudinal structur- 
ing on periods and damping times for kink modes in two- 
dimensional threads, and extracting conclusions about their im- 
plications for prominence seismology, our study aims at explain- 
ing the obtained results. For this reason, we have performed the 
analysis of the spatial distribution of perturbations and the en- 
ergy flux into the resonances. These two analyses, which are 
novel in the context of prominence oscillations, provide us with 
a comprehensive explanation for the obtained parametric results 
and add further insights to the physics of resonantly damped kink 
modes in two-dimensional equilibrium states. 

The paper is organised as follows. Section [2] describes 
the thread model, the relevant parameters introduced by its 
two-dimensional character, the linear MHD wave equations to 
be solved, and the numerical method used for that purpose. 
Section[3]presents our analysis and results. We first show com- 
putations of periods, damping times, and damping ratios as a 
function of the length of the thread and the density in the evacu- 
ated part of the tube. Implications for the determination of phys- 
ical parameters in prominences are discussed. Next, a qualita- 
tive explanation of the obtained results is given by analysing 
the spatial distribution of perturbations. Finally, we describe our 
energy analysis, that in combination with the spatial structure 
of eigenfunctions fully explains the obtained damping times. In 
Section|U our conclusions are presented. 

2. Thread model, linear MHD wave equations, and 
numerical method 

We consider an individual and isolated prominence thread in 
a gravity-free static equilibrium in the zero plasma-/? approxi- 
mation. The fine structure is modelled by means of a cylindri- 
cally symmetric flux tube of radius a and length L. In a sys- 
tem of cylindrical coordinates (r, ip, z) with the z-axis coinciding 
with the axis of the tube, the magnetic field is pointing in the 
z-direction and has a uniform field strength. Because of the as- 
sumed zero-/? approximation the density profile can be chosen 
arbitrarily. The non-uniform thread is then modelled as a den- 
sity enhancement with a two-dimensional distribution of den- 
sity, p(r,z) (see Fig. [TJ. The density distribution has two non- 
uniform layers, with length / and l z in the r- and z- directions, 
respectively. The first is introduced so as to study the resonant 
damping of oscillations. The second produces irrelevant phys- 
ical results, but enables us to avoid contact discontinuities and 
provides us with a continuous background density. Surface plots 
of the density distribution for different thread models are shown 
in Figure [2] where we have made use of the symmetry of the 
system in the r- and z- directions and only plot their positive 
values. The dense part of the tube with prominence conditions, 
i.e., the thread, has a density pf and occupies only part of the 
larger magnetic flux tube. It extends over a length Lthread in the 
z- direction. The rest of the internal part of the tube, with length 
L-Lthread in the longitudinal direction, is filled with plasma with 
a density p ev , with the subscript "ev" indicating the evacuated 
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Fig. 1. Schematic representation of the model configuration adopted in this work. The two-dimensional cylindrically symmetric 
line-tied structure consist of a dense part, the thread, with length L t h rea d and density pf, and two evacuated parts, with density p ev . 
Both regions are separated by a non-uniform layer of width l z in the longitudinal direction. Transverse non-uniformity is considered 
on a layer of width I. The structure is surrounded by plasma with coronal properties and density p c . 




Fig. 2. Density distribution in the (r, z)-plane in the domain r e [0, r max ], z e [0, L/2] for prominence thread models with non- 
uniform radial and longitudinal structuring. The threads are defined by a cool and dense part with density, pf and length Lthread 
and an evacuated part of the tube with density p ev and length L - L t h re ad- These two regions connect non-uniformly to the coronal 
surrounding medium with density p c . Depending on the value of p ev , three different situations are possible: (a) p ev = p c , (b) 
<p f , and (c)p ev <p c . 



part of the tube. All lengths a re normalised by taking a — 1 . In 
contrast to ISoler et ail (|2010i), the value of p ev can be different 
from the coronal density, p c , and may have values lower than p c 
up to a value equal to the filament density, pf , case in which the 
full tube is occupied with dense cool plasma and we recover the 
one-dimensional case. Both regions along the axis of the tube 
are connected by means of a non-uniform transitional layer of 
length l z to produce a smooth longitudinal profile. The density 
variation at this layer, with L t hread - h/2 < z < ^thread + can 
be expressed as 



A(Z) = Y 



Pf 



1 + ^|_h_^ 



Pf 



n(Z - Lthread) 



I, 



(1) 



for r < a - 1/2. As for the radial direction, the internal filament 
plasma, with density pf, is connected to the external coronal 
plasma, with density p c , by means of a non-uniform transitional 
layer of thickness /, defined in the interval [a - 1/2, a + 1/2], that 
can vary in between Ija = (homogeneous tube) and l/a — 2 
(fully non-unif orm tube). In contra st to the one-dimensional 
model used by lArregui et alj d2008l) the dense plasma of the 
thread does not occupy the full length of the tube, hence the 



radial variation of the plasma density for z < ^thread - and 
a - 1/2 < r < a + 1/2 is given by 



Pr(r) = — 



Pf 



p c \ . n(r-a) 
1 sin 



Pf 



/ 



(2) 



Radial non - unifor mity is not restricted to the dense part, as in 
Sole r et al. (2010), but can be present also in the evacuated part 
of the tube. As a consequence, there is an overlap region of radial 
and longitudinal non-uniform layers, for a - 1/2 < r < a + 1/2 
and Lthread - IJ2 < z < Lthi-ead + l z /2, where the density is given 
by 



Pr Z (r, Z) 



Pz(z) 



1 + 



Pz(z) 



Pz(z) 



sin ■ 



n(r — a) 



I 



(3) 



In the evacuated part of the tube, for z > £thread+'z/2 and a-l/2 < 
r < a + 1/2, the radial density profile is given by 



i \ P™ 
PM = — 



1 h 1 sin 



Pe 



Pf 



(4) 



Note that since the case p ev < p c is not excluded, the slope of 
the radial density profile in the evacuated part of the tube can be 
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positive (if p ev < p ev ), negative (if p ev > p ev ) or zero (if p ev = 
p ev )- Finally, for r > a + 1/2 we reach the coronal medium and 
p(r, z) — p c - The model adopted in this paper removes unrealistic 
discontinuities in the density distribution, considered in previous 
works, and allows the theoretical modelling of a large number of 
threads by simply considering different geometrical and physical 
parameter values for, e.g., the three different situations outlined 
in Figure [2] 

This paper is concerned with standing kink waves in promi- 
nence threads. To study small amplitude thread oscillations, we 
consider the linear resistive MHD wave equations for perturba- 
tions of the form f(r, z) ~ exp(i(a>t + mip)) with constant resistiv- 
ity, T]. Here m is the azimuthal wave-number and w = wr + iaty, 
the complex oscillatory frequency. For resonantly damped solu- 
tions, the real part of the frequency gives the period of the os- 
cillation, P - 2n/a>R, while the imaginary part is related to the 
damping time, tj = \/a>\. Magnetic diffusion is only included 
here to avoid the singularity of the MHD equations at the reso- 
nant position, but diffusion has no effect on the resonant damping 
time-scales. Oscillations are then governed by the following set 
of partial differential equations for the two components of the 
velocity perturbation, v r and v^, and the three components of the 
perturbed magnetic field, b r , b<p, and b z , 



iu)b r 



B_idbj_ _ dbz 

HP \ dz dr 
B (i m b db v 
pp \ r dz 
8 2 b r 
dr 2 



n dv r 
B ~z +T1 



nr , d b r 1 db r 



-rb r + 



dz 2 



r dr 



im b r 



itob, B = 



dvu, 



d 2 K 



dr 2 



irK + 



+ 2—b r - - 



. , B (dv r v r im 

iu)b 7 - -B I — — I 1 Vu 

dr r r 



+ 1 



d 2 ^ 1 db v 
dz 2 r dr 



d 2 b, m 2 
irb, 



dr 2 



d 2 b z 

dz 2 



1 db z 
r dr 



(5) 
(6) 

(7) 

(8) 

(9) 



Equations (J5J-®, together with the appropriate boundary con- 
ditions, define an eigenvalue problem for resonantly damped 
modes. As the plasma-/?=0, the slow mode is absent and there are 
no motions parallel to the equilibrium magnetic field, v z = 0. We 
further concentrate on perturbations with m = 1, which repre- 
sent kink waves that produce the transverse displacement of the 
axis of the tube. The MHD kink wave represents a wave mode 
with mixed fast and Alfven character, its Alfve nic nature being 
domi nant in and around the resonant position dGoossens et al.L 
2009). The problem can be further simplified by making use of 
the divergence-free condition for the perturbed magnetic field 



1 d(rb r ) 
r dr 



I dby db, 
r dip dz 



0, 



(10) 



which reduces the system of equations to be solved to four, upon 
expressing b^ in terms of b r and b z . Solutions to these equations 
are obtained by performing a normal mode analysis. Because 



of the complexity of the problem when a two-dimensional den- 
sity p(r,z) is considered, numerical solutions to the frequency 
and spatial structure of eigenfunctions in the (r,z)- plane are 
computed using PDE2D (Sewell, 2005), a general-purpose par- 
tial differential equation solver. The code uses finite elements 
and allows the use of non-uniformly distributed grids, which 
are needed to properly resolve the large gradients that arise in 
the vicinity of resonant positions. Different grid resolutions have 
been tested so as to assure the proper computation of the res- 
onant eigenfunctions. Magnetic dissipation has no influence on 
the damping of kink modes due to resonant absorption, a condi- 
tion that has to be checked in all numerical computations, but 
allows us to properly compute the spatial distribution of per- 
turbations in the resonance. We have made use of the symme- 
try of the system and solutions are computed in the domain 
r G [0, r max ], z e [0, L/2]. Non-uniform grids are used in both di- 
rections, to properly resolve the regions with r e [a - 1/2, a + l/2] 
and z e [ithread/2 - / z /2,L t h rea d/2 + l z /2]. The first is located so 
as to include the non-uniform transitional layer in the radial di- 
rection, while the second embraces the non-uniform transitional 
layer along the tube. As for the boundary conditions, we equate 
them to the spatial distribution of perturbations for the funda- 
mental kink mode. The two perturbed velocity components, v r 
and v,p and the compressive component of the perturbed mag- 
netic field, b z have vanishing longitudinal derivatives at z = 0, 
which corresponds to the apex of the flux tube, while they van- 
ish at z = L/2, because of the line-tying boundary condition at 
the photosphere. In the radial direction, they also have vanish- 
ing radial derivative at the axis of the tube, r = 0, while we 
impose the vanishing of the perturbed velocity components far 
away from the tube in the radial direction, hence (v r , v^) — > as 
r — > oo, a condition that is accomplished by setting the pertur- 
bations equal to zero at r — r max , where r max , the upper limit of 
the domain in the radial direction, has to be chosen to be suffi- 
ciently far to properly compute the drop-off rate of perturbations 
in the radial direction. In all our computations, we have consid- 
ered r max = 20a. 

3. Analysis and results 

Arre gui et al.1 (120081) . in their analysis of the damping of kink 
oscillations in one-dimensional thread models, showed that the 
parameters that determine the temporal attenuation of oscilla- 
tions are the density contrast, pf/p c and the width of the non- 
uniform transitional layer, I /a. The damping ratio is rather de- 
pendent on the first parameter for low values of the density 
contrast, but stops being dependent in the high contrast ratio 
regime, typical of prominence plasmas. The strongest influence 
comes from the width of the transitional layer, with the damp- 
ing time rapidly decreasing for increasing values of l/a. Before 
we deal with the additional parameters introduced by the lon- 
gitudinal density structuring in two-dimensional thread models, 
solutions to Equations ©-(O have been first obtained by us- 
ing a two-dimensional (2D) density distribution with Lthread = L. 
The purpose of these numerical experiments has been to check 
the correc t behaviour of the co de by reproducing the results ob- 
tained by lArregui et al.l (120081) . In addition, we have also con- 
sidered the magnetic Reynolds number, R m = v^aji], which 
should not affect the computed damping times, in the limit of 
large Reynolds numbers. Figure [3] displays the obtained results. 
The damping time of resonantly damped kink waves is indepen- 
dent of the magnetic Reynolds number, as long as this quantity is 
large enough for resonance absorption to be the operating damp- 
ing mechanism. This regime (see the plateau regions) is obtained 
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Fig. 3. Damping time, in units of the internal filament Alfven crossing time, tm = a/vAf, as function of the three relevant parameters 
for kink oscillations in one-dimensional thread models: (a) the magnetic Reynolds number, R m = Vf^a/rj, (b) the density contrast, 
Pf/p c , and (c) the transverse inhomogeneity length-scale, l/a. Solid lines indicate the ID solution, while symbols represent the 
numerical 2D solutions obtained with Lthread = L. In (a) and (c) the density contrast is pf /p c = 200. In (b), l/a - 0.2. In (b) and (c) 
R m = 10 6 . All computations have been performed in a two-dimensional grid with N r =401 and N z = 5 1 points, with 250 grid-points 
in the radial transitional layer. Lengths are normalised to a = 1 and L = 50a. 
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Fig. 4. (a)-(c): period, damping time, and damping ratio as a function of the length of the thread for thread models with l/a — 0.2 
and for two values of the density contrast, (d)-(f): the same quantities for thread models with pt/p c = 200 and for three values of the 
transverse inhomogeneity length-scale. In all figures a = 1, R m = 10 6 , p ev = Pc, and l z = a. Times are shown in units of the internal 
filament Alfven crossing time, T&f = a/vpj. Symbo ls correspond to fully numerical 2D computations while the different line styles 
represent the approximate 2D solutions obtained bv lSoler et al.l d2010t) . All computations have been performed in a two-dimensional 
grid with A^,- = 401 and = 51 points, with 250 grid-points in the resonant layer. Lengths are normalised to a = 1 and L = 50a. 



for different values of R m when different transitional layers are 
considered. Figure [3^ shows a perfect agreement between the 
ID results and the current computations using the 2D code. The 
perfect correspondence between ID and 2D computations for 
the damping time as a function of the density contrast in the ra- 
dial direction is shown in Figure fJJ). Finally, the most important 
parameter that determines the damping of transverse thread os- 
cillations is the width of the non-uniform transitional layer in the 
radial direction. The damping time strongly decreases when this 
parameter is increased, as can be seen in Figure[3j;, which again 
shows a very good agreement between the values computed in 
2D and the previous ID computations. These results were in 



agreement with previous works in the co ntext of the damping of 
coronal loop transverse oscillations (e.g. Goossens et all [ 1 992 ; 



iRuder man & Roberts! l200l iGoossens et al.U2002l) 



Once we are confident about the goodness of the code, we 
consider the new ingredients introduced by the two-dimensional 
nature of the prominence thread models considered in this work. 
These new ingredients are the length of the thread, Lthread, i.e., 
the length of the part of the magnetic flux tube filled with dense 
absorbing plasma, and the density in the evacuated part of the 
tube, p ev . The non-uniform transitional layer in density between 
both regions along the tube, l z , has an irrelevant effect on the os- 
cillatory period for low harmonics in the longitudinal direction, 
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as shown by iDfaz et al . (2008). Our comp utations (not shown 
here) confirm the finding by Diaz et ail (120081) on the negligible 



importance of this parameter for the oscillatory period and show 
a similar irrelevance concerning the damping time by resonant 
absorption. For this reason, we have concentrated our analysis 
on the remaining two parameters, L t h re ad and p ev . 



3.1. Periods and damping times 

We first consider the influence of the length of the thread on 
the period and damping of resonantly damped transverse thread 
oscillations. An in itial analysis on this subject was presented by 
ISoler et all d2010l) . who considered the thin tube and thin bound- 
ary approximations (TTTB) in a two-dimensional thread model 
with transverse inhomogeneity only in the dense part of the tube. 
Our analysis goes beyond the TTTB approximations by consid- 
ering a fully inhomogeneous two-dimensional density distribu- 
tion and combining its influence with that of the density in the 
evacuated part of the tube, p ev . 

We have first analysed the variation of period, damping time, 
and damping ratio, Td/P, by setting p ev = p e , so that we mimic 
the case studied by Sole r et all (l2010h analytically. We start with 
a fully filled tube and gradually decrease the length of the thread. 
The obtained results, for two values of the density contrast be- 
tween the filament and coronal plasma, are shown in Figures [4^- 
c. The period is strongly dependent on the length of the thread. It 
decreases by almost a factor of two when going from Lthread = L 
to Lthread = 0.1L. Figure [4^ also shows that the oscillatory pe- 
riod is almost independent of the density contrast, once this pa- 
rameter is large enough. A s for standing kink waves in one- 
dimensional thread models (Arre gui et all [2008). the kink fre- 
quency is a weighted mean of the internal and external Alfven 
frequencies. Regardless of the density contrast, the period is al- 
lowed to vary in a a narrow range determined by a factor that 
goes from V2 to 1, when going from pf/p c — 1 to pf/p c — > °°. 
For typical density contrasts in prominence plasmas, the pe- 
riod can be considered independent of the density contrast. The 
damping time produced by resonant absorption (Fig.|4p) also de- 
creases remarkably when the length of the cool and dense part 
of the tube is decreased. The decrease is also arou nd a factor 
of two in the considered range of values for Lthread- ISoler et aTl 
(2010) find that in the TTTB limit the dependence of the period 
and the damping time with the length of the thread is exactly the 
same, hence any influence on the damping ratio, Td/P, is can- 
celled out. Outside the TTTB approximations, we find that this 
is not the case (see Fig. |4j;), although the damping ratio is al- 
most independent of the length of thread and only for very short 
threads a slight increase in the damping ratio is found when fur- 
ther decreas ing this parameter . In Figure |4] we overplot results 
obtained by ISoler et ail d2010 ) by solving their dispersion rela- 
tion. We see that there is a very good agreement and the differ- 
ences, that are due to the simplifying assumptions of the analyt- 
ical treatment, are rather small. One can obs erve an anomalou s 
behaviour on the damping time computed bv lSoler et al.1 (12010). 
for small values of L t hread- This is due to the simplifying assump- 
tions considered to obtain the semi-analytic solution, that might 
not be entirely valid outside the long-wavelength limit. Overall, 
our results c onfirm the v alidity of the analytical approximations 
obtained bv lSoler et all (120161) concerning the influence of the 
length of the thread on periods and damping times. In physi- 
cal terms, the shortening of the length of the thread produces 
shorter period oscillations, since the physical system is equiv- 
alent to a fully filled tube, with the wavelength of oscillations 



replaced by a shorter effective wavelength. A physical explana- 
tion of the damping time dependence on the length of the thread 
is provided, by using energy arguments, in Sect. 13.31 

Similar conclusions can be extracted from the computa- 
tions we have performed by fixing the density contrast and for 
three different values of the width of the inhomogeneous layer. 
Figures |4ji-f show the obtained results. They clearly show how 
strongly the damping time and the damping ratio are influenced 
by the width of the transitional layer, also in 2D models, while 
the period of the oscillations is almost unaffected by the value 
of l/a. In view of the results displayed in Figure [4] we conclude 
that the length of the thread is a very important parameter. When 
allowing to vary from the limit of fully filled tube to 10% filled 
tube, periods and damping times are decreased as much as 50% 
percent. This is very relevant in connection to prominence seis- 
mology. We must note that, in principle, the length of the thread 
can be estimated directly from observations. However, the length 
of the supporting magnetic flux tube is much more difficult to es- 
timate, since its end points are usually unobservable. 

Our general density model enable us to analyse the effect of 
the density in the evacuated part of the tube on the oscillatory 
properties, as well. The study of th e influe nce of this parameter 
was not undertaken by Sole r et all (1201 Oh . and requires a fully 
numerical approach. We allow for p ev to be different from the 
coronal density, p c . When changing the density in the evacuated 
part of the tube, in addition to the radial non-uniform layer that 
connects the prominence material to the corona we are introduc- 
ing an additional radial non-uniform layer in between the evac- 
uated part of the tube and the corona. The density profile in this 
layer and its slope depend on the relative values of p ev and p c , 
as given in Equation and shown in Figs.|2p and c. Instead of 
analysing the effect of p ev on the oscillatory properties in a sepa- 
rate manner, we have selected three representative values for the 
length of the thread and have computed periods, damping times, 
and damping ratios as a function of the density in the evacuated 
part of the tube, measured in units of the coronal density. We 
have split our analysis into two parts. 

First, we consider that p c < p ev < pf, hence the parameter 
is allowed to vary in between the coronal and prominence den- 
sities. Fig.|5]displays the obtained results. For p ev = 200p c = pf, 
the tube is fully filled with cool and dense plasma. Then, as we 
decrease p ev , periods and damping times have a marked linear 
decrease. When 40% of the tube is filled with cool plasma, they 
decrease by about a 15%. When a 10% filled tube is consid- 
ered a decrease of up to a 50% is obtained. In physical terms, 
the period decrease when the density in the evacuated part of 
the tube is gradually decreased can be understood if we think 
about the different inertia of the system, with or without dense 
plasma at those locations along the tube. For explaining the dif- 
ferent damping time-scales, energy arguments that combine the 
energy of the mode and the energy flux into the resonance have 
to be considered, see Sect. 13. 31 The decrease in period and damp- 
ing time is very similar, but not exactly the same. For instance, 
Fig. [5}; shows that the damping ratio is slightly dependent on the 
density in the evacuated part of the tube, for all the considered 
values in the range p c < p ev < pf. 

Next, we have considered the possibility of the density in 
the evacuated part of the tube b eing lower than the coronal 
density. For instance, iDfaz et al (2002) considered a value of 
Pev = 0.6p c . The computations shown in Fig.|5]are extended to 
lower values of p ev . The obtained results are displayed in the in- 
set plots of Fig. |5]and show that the density in the evacuated part 
of the tube is irrelevant in relation to the period of the oscilla- 
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Fig. 5. Period, damping time, and damping ratio as a function of the density in the evacuated part of the tube for prominence threads 
with pf/pc = 200, a = 1, ^ m = 10 6 , l z = a, and l/a = 0.2 and three values of the length of the thread. The main plots are for 
Pc ^ Pev ^ Pi, while the inset plots correspond to the range 0.1p c < p ev < p c - Times are shown in units of the internal filament 
Alfven crossing time, t\{ = a/vAf- All computations have been performed in a two-dimensional grid with N r = 401 and N z = 51 
points, with 250 grid-points in the resonant layer. Lengths are normalised to a = 1 and L = 50a. 





Fig. 6. Longitudinal and radial dependence of the eigenfunctions transverse velocity component, v r , azimuthal velocity component, 
v,p, and compressive magnetic field component, b z , in prominence threads with pt /p c = 200, l/a - 0.2, a = 1, p ev = Pc, and 
R m = 10 6 , for different values of the length of the thread. All computations have been performed in a two-dimensional grid with 
N r = 401 and N z = 51 points, with 250 grid-points in the resonant layer. 



tions and the damping by resonant absorption in the considered 
range of values with 0.1p c < p ev < Pc- 

These results show that the density in the evacuated part of 
the tube is also a relevant parameter for prominence thread seis- 
mology, because of its effect on periods and damping times and 
the difficulty in being measured by direct means. When consid- 
ered in combination with effects due to the length of the thread, 
out computations enable us to perform a more accurate promi- 
nence seismology, applicable to a large number of threads. 

3.2. Spatial distribution of eigenfunctions 

Changes in the longitudinal density structuring have a direct 
impact on the oscillatory period, which is easy to understand, 
but also on the damping time by resonant absorption. Besides 



obtaining the parametric behaviour of kink mode periods and 
damping times as a function of the longitudinal density structur- 
ing, we aim to explain the obtained results. As a first step, we 
have analysed the spatial structure of eigenfunctions. The rel- 
evant perturbed quantities are the radial and azimuthal veloc- 
ity components, v r and v^, and the compressive component of 
the perturbed magnetic field, b z , directly related to the magnetic 
pressure perturbation, Pj = Bb z /fi. 

We first consider the influence of the length of the thread on 
the profiles of the eigenfunctions in the radial and longitudinal 
directions. Results are given in terms of the modulus of the com- 
plex eigenfunctions. Fig. [6] shows one-dimensional cuts along 
the longitudinal and radial directions of the eigenfunctions for 
different values of the length of the thread. The longitudinal pro- 
files are shown at the axis (r = 0) for v r , and at the mean radius 
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of the tube (r = a) for v v and b z . In the longitudinal direction all 
three eigenfunctions display a trigonometric dependence with z, 
when Lthread = L, i.e., the case that mimics the one-dimensional 
thread. When the length of the thread is decreased several inter- 
esting effects occur. First, both perturbed velocity components 
display a slightly improved confinement. The maximum values 
of the velocity perturbations still occur in the dense part of the 
tube, but the drop-off rate changes, becoming almost linear out- 
side the thread so that they satisfy the boundary condition at the 
foot-point of the tube. When eigenfunctions are normalised to 
the value of |v r ] at the apex of the tube, it is seen that the decrease 
of the length of the thread produces larger amplitudes of the az- 
imuthal velocity component (see Fig.|6j)), related to the Alfvenic 
character of the mode. This amplitude almost doubles its value 
when going from L t hread = L to L t hread = 0.04L. The z-component 
of the magnetic field perturbation gives an indication of the com- 
pressibility of the normal mode. Our results indicate that when 
the length of the thread is decreased, the longitudinal profile of b z 
(Fig-ISt) becomes strictly confined to the dense part of the tube, 
where it reaches its maximum value. As the thread gets shorter, 
the maximum value of b z increases, hence the compressibility 
becomes larger at the apex of the tube for shorter threads oscil- 
lations, though the kink mode remains being an almost incom- 
pressible wave mode. Kink modes in fully filled magnetic flux 
tubes are almost incompressible. Longitudinal density structur- 
ing by the inclusion of a dense central part increases (decreases) 
kink mode compressibility in the dense (evacuated) parts of 
the tube, in comparison to the one-dimensional flux tube kink 
modes. 

A close look at the spatial structure of eigenfunctions in the 
radial direction (Figs. |6ji-f) provides us with a qualitative expla- 
nation on why the change of the length of the thread affects the 
damping times co mputed in the pr evious subsection (and also 
those presented bv lSoler et al.ll2.Q10» . First, the change in the pe- 
riod produced by the change in the length of the thread affects 
the resonant position in the radial direction and hence the damp- 
ing time. However, the slope of the density profile in the tran- 
sitional layer is very large, because of the high density contrast 
ratio typical of filament threads and this effect is not very impor- 
tant. Second, the eigenfunctions in the resonant layer are affected 
by the value of the length of the thread. Both effects are clearly 
seen in Figs. [6ji-f, where the radial dependence of v r , v v , and b z 
is plotted for different thread lengths. We first note that the per- 
turbed velocity components at the resonance have a larger am- 
plitude the shorter the length of the thread. We interpret this as 
an indication of an improved efficiency of the resonant damping 
mechanism. This is particularly clear if we look at the azimuthal 
velocity profile in the resonant layer (Fig. |6^). Now, not only 
the increase of the amplitude for shorter lengths is evident, but 
also the transverse length scale is seen to decrease when shorter 
threads are considered, hence thinner resonance widths are ob- 
tained. This detail of the resonance also shows a slight shift to- 
wards the left hand side which is due to the change of the reso- 
nant position for different lengths of the threads (hence different 
oscillatory periods). Although the resonances are not exactly lo- 
cated at ta = a, our numerical res ults indicate that this approxi- 
mation, used bv lSoler et alj d2010l) . is fully justified. Finally, the 
radial profile of the compressive component of the magnetic field 
perturbation (Fig. |6f) shows the increase in its amplitude at the 
resonant position mentioned above, when considering shorter 
threads. This is another indication of the improved efficiency of 
resonant damping. Notice that although our study is limited to 
linear MHD waves, in a realistic situation perturbations inside 
the resonant layer become nonlinear, as shown in studies by e.g., 



iTerradas et ail (l2008ah : Iciack et alJ J2009h : iBallai & Rudermanl 
(1201 lh . 

We have next examined the spatial distribution of the rel- 
evant perturbed quantities as a function of the density in the 
evacuated part of the tube. Fig. [7] shows the obtained profiles 
in the longitudinal and radial directions for a fixed value for 
the length of the thread. As before, the longitudinal profiles are 
shown at the axis (r = 0) for v r , and at the mean radius of the 
tube (r - a) for and b z . Our results indicate that the den- 
sity in the evacuated part of the tube has also a direct impact 
on the radial and longitudinal profiles of eigenfunctions. For the 
longitudinal profiles, as we decrease the value of p ev from pf, 
we have a slightly improved confinement of the velocity pertur- 
bations and a strict confinement of the longitudinal component 
of the perturbed magnetic field to the dense part of the tube (see 
Figs.|7^-c). The amplitude of b z at the apex of the tube increases, 
while in the evacuated part of the tube, and for a fixed length 
of the thread, compressibility decreases as p ev is decreased (see 
Fig. [7}:). Taking a look at the eigenfunctions in the radial di- 
rection (Figs. |7ji-f), the amplitude of both perturbed velocity 
components increases when the density in the evacuated part of 
the tube is decreased. The increase in the resonance peak and 
the shortening of the transverse spatial scale are mainly evident 
for the azimuthal velocity component, which gives its resonantly 
damped and Alfvenic character to the mode. We find an increase 
in the compressibility of the normal mode (Fig.[7]f) in the dense 
prominence plasma region as the density in the evacuated part of 
the tube is decreased. 

Overall, the decrease of the density in the evacuated part 
of the tube, starting from a fully filled tube, produces similar 
qualitative effects on the radial and longitudinal profiles for the 
eigenfunctions as the ones produced by the shortening of the 
length of the thread. In the parameter range studied in this work 
those effects are quantitatively more important in the case of the 
changes of the length of the thread. The properties of the spa- 
tial structure of eigenfunctions give us a qualitative explanation 
on why changes in the longitudinal density structuring in promi- 
nence threads have a significant effect on the damping time of 
kink oscillations in two-dimensional thread models. 

3.3. Energy analysis 

Our analysis of the spatial structure of eigenfunctions indicates 
that the decrease of both the length of the thread and the den- 
sity in the evacuated part of the tube produce a strengthening of 
the resonance absorption process that becomes apparent through 
the appearance of more pronounced resonant profiles and thin- 
ner resonance widths in the perturbed velocity components at 
the resonance. This would explain the marked decrease of the 
damping times found in Section lXTl when varying these two pa- 
rameters. This result shows that resonant damping may strongly 
depend on the details of the longitudinal density structuring in 
rather general two-dimensional density models, such as the ones 
used here to model prominence threads. 

In order to obtain a quantitative explanation an energy anal- 
ysis is carried out. Our analysis involves the energy of the kink 
mode and the energy flux into the resonance. The energy of the 
mode can directly be computed from the numerical eigenfunc- 
tions as 

E= l -(pv 2 +b 2 ). (11) 

This expression has to be integrated over the entire (r, z)- 
plane except for the resonance layer, easily recognisable from 
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Fig. 7. Longitudinal and radial dependence of the eigenfunctions transverse velocity component, v r , azimuthal velocity component, 
v^, and compressive magnetic field component, b z in prominence threads with pf/p c = 200, l/a = 0.2, a = 1, L t hread = 0.2L, 
and R m = 10 6 for different values of the density in the evacuated part of the tube. All computations have been performed in a 
two-dimensional grid with N r = 401 and N z = 51 points, with 250 grid-points in the resonant layer. 



the eigenfunctions displayed in Figs. [6] and [7] as this region 
would include a contribution from the Alfven waves. In one- 
dimensional equilibrium models, the energy flux into the res- 
onance i s proportiona l to the magnetic pressure pertur bation 
squared ( Andr ies"et al.L 120001; lAndries & Goossensl 12001). a re- 
sult that was used by Arregui et al.l (l2007bl) to analyse the influ- 
ence of the internal structuring of coronal loops on the damping 
by comparing the efficiency of the process at internal and ex- 
ternal layers. In two-dimensional equilibrium states, the energy 
flux absorbed at a particular field line is proportional to the over- 
lap integral between Pj{rA,z), the profile of Pj along the tube 
at the resonant position, and t he resonant Alfven eigen functions 
dThompson & Wright! 119931: iTirrv & Goossensl Il995l) . This is 
why the longitudinal profiles of b z and v v , and their modification 
due to changes in the longitudinal density distribution are so rel- 
evant in determining kink mode damping times. A general math- 
emat ical expression, valid in the three-dimensional case, is given 
bv lWright & Thompsonl (Il994l) . When adapted to our cylindri- 
cal equilibrium, this expression gives the energy flux at the res- 
onance per unit <p in the form 



F = r A J [S r (r A )-S r (rX)]dz 

nm 2 B 2 ( f(f> b z dzf (da> A 

— ^ ^ I 

4|tg r A ffipfidz \ dr 



(12) 



In this expression, S r (r A ) and S r (r A ) represent the values of the 
radial component of the time-averaged Poynting vector to the 
left and to the right of the resonance position, that we approx- 
imate by r A — a, <p is the Alfven eigenfunction along the tube 
at the resonant position, and cop,{r) the Alfven continuum fre- 
quency. In order to evaluate the integrals as well as the slope of 



the Alfven continuum at the resonant position in expression (fT2l . 
we have first solved the following equation for the Alfven con- 
tinuum modes 



(r A ,z) u A (r A ) 



dz 2 



v l A {r A ,z) 



</>(r A ,Z) = Q, 



(13) 



with boundary conditions 

^!( z = O) = O,0(z = L/2) = O. 
dz 

For each value of this equation is solved for different val- 
ues for the relevant parameters Lthread and p ev , that in turn define 
different profiles for v 2 A (r A ,z). The Alfven continua arise from 
repeating the procedure for different values of r A in the range 
[a - 1/2, a + 1/2]. Once this is done, the slope of the Alfven con- 
tinua at the resonant position and the Alfven eigenfunctions that 
correspond to each thread model can be evaluated. By using the 
numerically computed profiles for b z (r^,z), in order to evaluate 
the overlap integral in Equation 111 2b . and multiplying by 27r, the 
required total energy flux is obtained. The ratio of the energy 
of the kink mode to the time-averaged energy flux into the res- 
onance gives the required damping time, for every considered 
value of Lthread and pev 

Figure [8] shows the obtained results. They have been nor- 
malised to the value of the damping time for the fully filled 
tube. We see an excellent agreement between the damping times 
computed through the energy analysis explained in detail above 
and the numerically computed ones for both cases in which we 
change the length of the thread and the density in the evacu- 
ated part of the tube. The results obtained in Section 13.11 can 
therefore be explained in terms of the energetics of the modes 
and the resonant energy transfer. These results also show the 



10 



I. Arregui et al.: Magnetohydrodynamic kink waves in two-dimensional non-uniform prominence threads 



(a) 



1.2 
1.0 

0.8 

0.6 

0.4 

0.2 




(a) 



0.2 



0.4 



0.6 



0.8 



1.0 




00 



1.2 
1.0 
0.8 
0.6 

0.4 
0.2 







50 



100 



150 



200 



(b) 




Fig. 8. (a) Damping time, normalised to the full tube damping 
time, as a function of the length of the thread, for p ev = p c . 
(b) Damping time, normalised to the full tube damping time, as 
a function of the density in the evacuated part of the tube, for 
^thread = 0.2L. In both figures lines correspond to the numeri- 
cally computed damping times and the symbols are the values 
obtained through the energy analysis described in Section [33l 



accuracy and utility of the analytical expression derived by 
Wright & Thompson (1 1 9941) for the time-averaged energy flux 
in terms of eigenfunctions. Furthermore, we have just demon- 
strated that damping time estimates can be obtained, using en- 
ergy arguments, without the need to solve the full non-uniform 
problem, but instead by solving the much more simpler piece- 
wise uniform problem in order to obtain the real part of the fre- 
quency and the longitudinal profiles for the eigenfunctions for 
kink modes in combination with the non-uniform computations 
for the Alfven continuum modes, using this information in com- 
bination with Equations (TTTI) and 111 2b . This is possible because 
the real part of the frequency and the longitudinal profiles of the 
eigenfunctions are only slightly affected by resonant coupling, 
and because the energy in Equation (TTTI) has to be integrated ex- 
cluding the resonant layer. Of course, the full numerical solution 
provides us with more accurate results, but as Fig. [8] illustrates 
the agreement between both methods is excellent. 

A convenient way to further understand the energy flux at the 
resonance and the dynamics of resonantly damped kink modes is 
to compute the spatial distribution of the time-averaged Poynting 
flux in our two-dimensional domain by making use of the gen- 
eral expression 



< S >= -Re(E x b*), 



(14) 



Fig. 9. (a) Spatial distribution in the (r, z)-plane of the time- 
averaged Poynting flux given by Equations ( TT3T > and (fT~6b . (b) 
Ohmic heating distribution, in arbitrary units. These results cor- 
respond to the transverse oscillation of a prominence thread with 
Pf/Pc - 200, pev = Pa 1 1 a = 0.2, a = 1, and Lthread = 0.2L. The 
magnetic Reynolds number is R m — 10 6 . 

with E = -(vxB) + ?7j the perturbed electric field, b*=(b* r , b* b* z ) 
the complex conjugate of the perturbed magnetic field, and j = 
(V x b)/p the current. In contrast to the energy analysis above, 
we now make use of the full numerical computations. In terms 
of the wave fields analysed in Section [3721 the two components 
of the time-averaged Poynting vector in the (r, z)- plane can be 
cast as 



<S r > 



<s z > = 



--Re(v r fop + -^Re 

2 /j ~ 2/j. 1 



dz dr 



r 



-b r 



(15) 



dby 
~dz 



ldb r 


db z \ 


\dz 


dr J 



(16) 



By making use of the divergence-free condition given by 
Equation ( ITOb to compute b^, we have produced an example of 
the two-dimensional distribution of the time-averaged Poynting 
flux around the resonant layer, for a partially filled thread. 

The arrow plot in Fig. [9^ shows that energy is fed into the 
resonant layer by concentrating it over the dense thread sec- 
tion. The amount of energy flux in t he radial direction is de- 
termined by the jump in < S T > (see IStenuit et all 1 1999L for a 
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one-dimensional example in ideal MHD). The resistive layer is 
on a scale where both the ideal and resistive terms of the per- 
turbed electric field can be of similar magnitude, but each of 
them could be dominant over different regions of the domain. 
For instance, the jump in < S r > is determined by the ideal 
term in Equation dT31 ) which is dominant in the thread region 
and zero in the evacuated part of the tube. Hence the energy in- 
flow towards the resonance in the radial direction comes from 
the interior of the tube and is determined by the spatial profiles 
of the transverse velocity component and the compressive mag- 
netic field perturbation. Note that the amplitude of the kink mode 
is smaller outside the tube compared to inside and also that the 
evacuated part of the tube is much less dense than the thread. 

Once in the layer, the energy flow is diverted along the field 
lines in a manner determined by < S z >, with a small radial con- 
tribution that is due to the resistive terms in Equation tl5i . Such 
as corresponds to the Alfvenic character of the mode inside the 
resonant layer, the dominant term in Equation ( TToT l involves the 
azimuthal velocity and magnetic field perturbations, v^Z?* . This 
quantity decreases linearly along the field lines, producing the 
lessening in the parallel Poynting flux towards the foot-point. 
In our example, the radial non-uniform layer is restricted to the 
dense part of the flux tube, since p ev = p c . However, energy con- 
centrated in the resonant layer of the thread because of resonant 
wave damping can flow along the field lines and, eventually, sup- 
ply heating in the evacuated region, where field aligned currents 
are dominant. Although the energy inflow into the resonance, 
given by < S r >, and its subsequent divergence along the field 
lines, given by < S z >, are mostly determined by ideal terms 
in Equations $15[ and ( TToT ). this does not mean that resistivity is 
not important. For instance, the amount of heating, in the form 
of Ohmic dissipation, will be determined by those currents, re- 
sistivity, and their spatial distribution. For this particular case, 
heating is distributed in a constant manner in the evacuated part 
of the tube, even if there is no resonant layer in that region (see 
Fig. |9h). 

4. Summary and Conclusions 

Quiescent filament fine structures are only partially filled with 
cold and dense absorbing material. The length of the threads 
can in principle be measured in events showing transverse os- 
cillations, provided the lifetime of threads is sufficiently large 
compared to the oscillatory period. The length of the supporting 
magnetic flux tubes are however much larger, and cannot be ob- 
served. Density measurements, both in the thread as in the evac- 
uated part of the supporting magnetic tube, are also challenging 
from the observational point of view. It is therefore important 
to quantify the variations on wave properties due to changes in 
these equilibrium parameters if we aim to perform an accurate 
prominence seismology. It is essential to have computations of 
periods and damping times for a wide range of thread models 
to include regimes in which the applicability of simple analyti- 
cal models could be of limited extent. For this reason we have 
computed the oscillatory properties of resonantly damped trans- 
verse kink oscillations in rather general two-dimensional fully 
non-uniform prominence thread models. This allows for a broad 
range of prominence threads with very different physical con- 
ditions to be modelled and their oscillatory properties charac- 
terised. 

The length of the thread and the density in the evacuated part 
of the tube define their longitudinal density structuring. We find 
that the length of the thread strongly influences the period and 
damping time of transverse kink oscillations, while the damping 



ratio is rather insensitive to this parameter. These resu lts confirm 
the val idity of the analytical approximations made bv lSoler et alJ 
d2010h . In addition, our modelling has allowed us to identify 
a new physical parameter with seismological implications, the 
density in the evacuated part of the thread. This quantity also in- 
fluences periods and damping times, and to a lesser extent damp- 
ing ratios, and must be taken into account in the inversion of 
physical parameters in the context of prominence seismology. 

Currently available inversion scheme s for one-dim e nsiona l 
coronal loops a nd prominence threads (|Arregui et all 12 007 a ; 
iGoossens et a I2008t lArregui & Ballesteii 1201 Cj ISoler et all 
2010) make use of observed periods and damping ratios. The 
first, influence the inferred values for the Alfven speed, while 
the second determine the transverse density structuring. Based 
on our results, we can conclude that ignorance on the length of 
the thread, the length of the supporting magnetic flux tube, and 
the density in the evacuated part of the tube will have a signif- 
icant impact on the inferred values for the Alfven speed (hence 
magnetic field strength) in the thread, depending on whether we 
use those one-dimensional inversion schemes or the results from 
two-dimensional models here obtained. On the contrary, because 
of the smaller sensitivity of the damping ratio to changes in the 
longitudinal density structuring, seismological estimates of the 
transverse density structuring will be less affected by our igno- 
rance about the longitudinal density structuring of prominence 
threads. 

Our study provides additional insight to the physics of res- 
onantly damped kink modes in two-dimension al equilibrium 
states, by extending prev ious applications (e.g, lAndries et all 
2005: lArregui et all l2005h to more complex non-separable den- 
sity distributions. It also provides an example of the methods 
and uses of combining the information from the spatial distribu- 
tion of eigenfunctions with that obtained from energy arguments. 
In particular, our energy analysis has allowed us to explain the 
decrease in damping times for shorter thread lengths found by 
ISoler et ail (l2010h . The length of the thread influences the en- 
ergy of the kink mode, and hence its oscillatory period, but also 
affects the damping by resonant absorption, through the energy 
flux into the resonance. In an analogous way, the value of the 
density in the evacuated part of the tube also determines periods 
and damping times, since both the energy of the kink mode and 
the energy flux into the resonance vary. This means that changes 
in the equilibrium configuration in a non-resonant direction pro- 
duce variations in the damping properties of kink modes, a result 
that was qualitatively explained by a detailed examination of the 
radial and longitudinal profiles of the eigenfunctions. Both the 
shortening of the length of the thread and the decrease of the 
density in the evacuated part of the tube produce more marked 
resonances, with the amplitude of the velocity perturbations at 
the resonance and the compressibility of the mode in the thread 
being larger. Inside the resonant layers shorter transverse spa- 
tial scales for the Alfvenic velocity component are obtained. In 
combination with the analysis of the energy of the kink modes 
and the energy flux into the resonances a quantitative explana- 
tion was obtained for both the dampin g properties obtained in 
our study and those in lSoler et al.l d2010h . 

The damping of kink oscillations in two-dimensional fully 
non-uniform equilibrium configurations can be computed by us- 
ing energy arguments together with the solution of simpler prob- 
lems for kink mode and Alfven continuum modes. This aspect 
is worth to be considered in future studies of resonant absorp- 
tion in 2D/3D models of solar atmospheric magnetic structures 
involving changes of equilibrium parameters that affect the den- 
sity structuring in a non-resonant direction. The use of energy 
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arguments would allow to have a first indication about how im- 
portant a given parameter that modifies the equilibrium in a non- 
resonant direction is, while avoiding to solve the full problem 
until we are interested in the details. 

The two-dimensional distribution of the time-averaged 
Poynting flux shows how energy is fed into the resonance and 
subsequently flows along the field lines by properties of wave 
fields associated to ideal processes. This is the reason why reso- 
nant damping is a mechanism for wave energy transfer in which 
time-scales are independent of resistivity, in the limit of high 
magnetic Reynolds numbers. Magnetic diffusion plays its role 
once energy is concentrated at small spatial scales, by provid- 
ing heating at locations that, as in our example for a partially 
filled thread, are distributed in regions where no resonant layer 
is present. This result offers additional insights to the dynamics 
of resonantly damped kink modes and the heating of prominence 
plasmas by wave transformation processes. It must be consid- 
ered in detail and extended to density models relevant to other 
solar atmospheric structures. 
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